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ON IDENTIFIABILITY OF FLEXIBLE STRUCTURE 
PARAMETERS 

By 

S. M. Joshi^ a^d G. L. Goglia^ 

ABSTRACT 

This report investigates the identifiability of modal parameters of 
flexible structures. Expressions are derived for Cramer-Rao bounds 

for the modal parameters, that is, frequencies, damping ratios and mode 
shapes or slopes. The optimal initial state, which maximizes the trace of 
the Fisher information matrix in the absence of persistent input, is obtain- 
ed. The concepts discussed are applied to a finite-element model of the 122 
meter hoop/column antenna. The numerical results show that the identifia- 
bility of the structural frequencies is excellent, followed by that of the 
danping ratios and the mode-slopes. 

INTRODUCTION 

The design of a controller for large flexible space structures requires 
reasonable knowledge of the structural mode parameters. Since it is usually 
difficult to obtain accurate knowledge of the parameters by apriori model- 
ing, it would be necessary to estimate the parameters in orbit. The iden- 
tifiability of the parameters is investigated in this report using the Cra- 
mer-Rao lower bound. The expressions for Fisher information matrix and the 
optimal initial state for identification are obtained in the next section, 
followed by numerical results of the application of these concepts to the 
122m hoop/ column antenna. 

^Research Associate Professor, Department of Mechanical Engineering and 
Mechanics, Old Dominion University, Norfolk, Virginia 23508 

^ Eminent Professor and Chairman, Department of Mechanical Engineering and 
Mechanics, Old Dominion University, Norfolk, Virginia 23508 
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IDENTIFIABILITY OF STRUCTURAL PARAMETERS 

The objective of this report is to investigate the identif iability of 
the parameters, i.e., mo<lal frequencies, damping ratios and mode shapes (or 
slopes) at the sensor locations. Since most established parameter estima~ 
tion methods use discretized sensor outputs, consider the discrete- time, 
constant-coefficient system: 


x(k->-l) 

“ Fx(k) + Gu(k) + v(k) 

(1) 

y(k) - 

Cx(k) + Du(k) 

(2) 

z(k) ■ 

y(k) + w(k) 

(3) 


where x, u, y represent nxl, mxl, £xl state, input and output vectors, v 
and w represent zero-mean, mutually uncorrelated white noise processes, and 
4 is the observation (sensor output vector). F, G, C, D are appropriately 
dimensioned matrices. If x(0), v(k), w(k) are Gaussian, the conditional 
density f(Z^/6) is Gaussian, 

where 

Zn “ {z(0),...,z(N)} (4) 

and 6 is the pxl vector of parameters to be estimated. Assuming v ■ 0, 

ff(Zjj/0) ■ Cl exp - j^l/2 {z(k) - y(k) ^ w”^ z(k) - y(k)[] (5) 

where Ci is a constant and W ■ cov. {w(k)l. (If the apriori density 
of 6, f(6), is Gaussian, the unconditional density f(Z^) will also be 
Gaussian) . 

If 9 is any absolutely unbiased estimate of 9 based on Z^, it can 
be proved that (ref. 1) 

E [(9 -9)(9 -9)’^/9] 


where 


E j LL- f (Zj,/9)/9 
( 39 ^ 


( 7 ) 


2 


It can be proved from equations (5 and 7) that 

^ w"^ ^y(h) 

39 36 


fel!- la 

"OOff QO/lUTi; 

( 8 ) 


J is the Fisher information matrix^ and J~^ is the Cramer-Rao lower bound 
on the covariance of the parameter estimation error. Thus ith diagonal ele- 
ment 1 b the lower bound on the variance of error in estimating 6^ , 

It is obvious from equation (8) that, J tends to equal a constant value 
as N-»- 00 , in the absence of persistent inputs to the system (provided that 
the system is asymptotically stable). J is nonzero only if the initial 
state is nonzero. If persistent inputs are present, J 4 - oe jb N oo, since 
the information available keeps growing. If all the parameters are 
identifiable, the information matrix is nonsingular. An understanding of 
the degree of identifiability of different parameters can be obtained by 
comparing the Cramer-Rao stavidard deviations for those parameters, obtained 
using nonzero initial state and zero inputs. 

Considerable work has been done on designing optimal inputs for identi- 
fication. However, for investigating the relative degree of identifiability 
of various parameters, zero input can be assumed. If nonzero inputs are 
used during identification, accurate measurements of the inputs are requir- 
ed. In addition, to be of any consequence, the inputs will have to have 
reasonably large magnitudes. This could be particularly objectionable for a 
LSS, as large magnitude inputs may cause excessive elastic motion which 
would jeopardize the structural integrity of LSS. Also, nonzero input (e.g. 
torque in the case of the antenna) would require the estimation of the mode- 
slope at the location of the actuator, which may add to the number of para- 
meters to be estimated. For these reasons, zero inputs are used in this 
study. The identification procedure could perhaps proceed as follows: af- 

ter the deployment of the antenna, one or more torquers could be used to 
create nonzero modal amplitudes. The torquers would then be turned off and 
the sensor data would be collected for N samples. The data could be used 
(on-line or off-line) in an identification algorithm to obtain the 
estimates. Control system gains could subsequently be computed, and the 
control system turned on. 
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OPTIMAL INITIAL STATE FOR IDENTIFICATION 


The sensitivity of y(k) to the ith parameter, is given by ^assutr- 
ing u ■ 0) 


3y(k) _ 9C 

afl aTT x(k) + C 


ae 


aT. 

X 


3x(k) 

a^; 

X 


(9) 


The sensitivity state vector is given by 


axXk+i) _ aF „ ax(k) 

Sr. 5T, ^ 55T“ 


( 10 ) 


The output sensitivity matrix is 


ay(k) 

ae 


3 yf kl ay(k) . ...» ay(k) 
aOj, 302 ae 


^Ej x^(k),...,EpX^(k) j 


J^xp 


wtiere the (p-t-1) n ^ dimensional, augmented state vector 


x^*^(k) -Ix^Ck), 3x(k) 3x(k) ^ 

L 99, 99 

1 P 


] 


(12) 


and are appropriate coefficient matrices. For zero input, the 

augmen^.ed state vector can be written as 


X (k) ■ * (k,0) X (0) 
a a 


(13) 


where ^(k,0) is the transition matrix for the augmented scate vector. 


Since 


x/(0) 


[x\0), o’^] 


( 14 ) 
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Equation (13) can be written aa 

X (k) - y(k) x(0) (15) 

a 

where Y is the appropriate submatrix of 4. Therefore, 

3y(k) ■|Ei(k) x(0), E 9 (k) x(0) E (k) x(0) L (16) 

ae p **p 

where E . (k) - E .t (k) 

1 1 

Substituting equation (16) into equation (8), the (i,j) th element of J is 
given by: 


* I " A T -1 A I 

J.. - x(0)< E E/(k) W E.(k)>x(0) (17) 

^ (k-0 ^ ) 

One definition of optimal initial state is the x(0) which maximizes 
the trace of J when constrained to be within a hyperellipsoid given by 

T 

X (0) S x(0) < d (18) 

T 

where S » S >0 and d > 0, 

Theorem 1. - The x(0) which maximizes Tr(j) subject to (18) is the 

eigenvector of S~^H corresponding to its maximum eigenvalue, which lies 
on the surface of the hyperellipsoid in (18), «4iere 

H “ " J 
i-1 

■i 

Proof . - The proof of this standard optimization problem is well known. 
Optimal Initial State for LSS. - The structural equations of motion for a 
flexible LSS are uncoupled. The equation of motion (assuning zero input) 
for the ith mode (i ■ 1, 2,.., n^) is: 


q*i + 2 p . w. q. + oj? q.- 0 


( 20 ) 
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The Jt X 1 output vector due to elastic notion is given by (asstsning rate 
gyros are used) 


y ■ jl , ♦ i’i 


( 21 ) 


The parameters to be estimated for the ith mode are: 


®i "If i» ‘"i* ^li» *2i»*’*»*Ai] 


( 22 ) 


Thus there are (t-f2) parameters per mode, or a total of n^(£-^2) 
parameters. Since the parameters appear explicitly in the continuous equa* 
tions of motion, it is more convenient to first write the continuous equa- 
tions for the sensitivity states, and then to discretize them. The sensi- 
tivity state equations for the ith mode with respect to p is: 


i + 2P . / ^^i \ + ^ *^i\ + 2m . q. ■ 0 


8p 


(23) 


(A similar equation can be written with respect to m^). In the absence 
of persistent input, the sensitivity with respect to the mode slopes is 
zero. Denoting 9^^ ■ [p^^, ^1 output sensitivity matrix 

or the ith mode is 


" r LZ- iz iz* • • • » iz 1 

39^ |^3p. 3(1) 3(|)]. j^ * ^*■^2) 


(24) 


t 

The output sensitivity terms on the right-hand side of equation (24) can 
be easily obtained from equation (21) and equation (24) can be rewritten 


as: 


3y fE-.x ,E.-X E. . _x 1 

g0 * I *i i2 a^ , . . . , i,A +2 a£ J 


where 


[ • 3q. 3q. 3q. 8q-l 

*11* Ji* Jli* ^ 

^ ^ 3p, 3p. 3m. SoTT 

1 1 1 iJ 


(2 5) 


(2 6 ) 
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using equation (26) in equation (8) and simplifying, the trace of the (i, 
i)th submatrix of J (i**!, 1, . . ,nq) can be written in the following form 
(since the initial conditions corresponding to the senriitivlty state varl~ 
ables are zero): 

Trj^ ■ (0) H xj (0) (2 7) 

where H. depends upon the state transition .matrix for equation (26) and W, 

^ • T 

and x^ ■ (q^q*) • Therefore, the initial state x^(0) for the ith 

mode, which will maximize Tr^ subject to the con.^itraint 

x^'^(O) S^x^(O) < (28) 

T 

(where S. ■ 0, d. > o) is the eigenvector corresponding to the maxi* 

mum eigenvalue of with a norm such that equality is satisfied in 

equation 28. 

The computation of the optimal initial state requires the knowledge of 
the parameters. The apriori (approximate) information about the parameters 
could be used for this purpose. Also, it will be necessary to compute the 
inputs in order to achieve the optimal initial state. This could be 
approximately done using the controllability matrix based on apriori para- 
meter information. 


NUMERICAL RESULTS 

In order to investigate the identifiability of various modal para- 
meters, Cramer-Rao (C-U) bounds were computed for the hoop-column an- 
tenna. This was done by first creating nonzero initial state using the 
torquer at actuator location no. 1, and then propagating the sensitivity 
state vectors. Two, 3 axis rate gyros were assumed at locations 1 and 2 
(fig. 2). Only the first 10 modes were considered (See Table l) . The para- 
meters to be estimated for the ith mode are: p., (n., 4*x^., 4>y,., 4>y„., 

1 i li 2i li 2i 
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f 


(P 0 ^ * since the iOth node has a natural frequency of 9.78 rad/sec 

*1 i 

1.5 Hz), the aampling frequency waa choaen to be 4 Hz (or sampling period 
- 0.25 see.). Since the (P<« )"^ for the slowest (first) mode is approxi- 
mately 130 sec., the data length for the compi.?tation was taken to be that 
value. The sensor noise in the rate gyrou asstined to have a standard 
deviation intensity (continuous- time) of 1 arc-sec/sec. (The C-R bounds for 
other noise levels can be obtained simply by multiplying by the appropriate 
factor) . 

An examination of the magnitudes of the mode- slopes for each mode cor- 
responding to X, y, z axes revealed the contributions of each mode to th e 
gyro outputs. These contributions were divided into three groups: strong 
(corresponding to large magnitudes of mode-slopet!i) , median, and weak 
(small magnitudes of mode-slopes). Table II shows this classification, 
where asterisk (*), tic-marc (/) and blank are used to denote strong, medi - 
um and weak contributions. It is sensible to try to identify only those 
modes using a particular axis gyro output, which make either strong or med- 
ium contribution along that gyro axis. Thus x-axis and z-axis gyro outputs 
are used to identify modes 1, 3, 4, 5, 6 , 8 , 10, while y-axis gyro outputs 
are used to identify modes 2, 7, and 9. 

Nonzero initial states were attained using torquers (along appropriate 
axis) at actuator location no. 1. The torquers were assumed to generate 
rectangular- waves with the same frequency as the lowest- frequency mode to be 
estimated (for that axis), with nominal amplitude of 10 ft-lbs. (optimal 
initial conditions were not used). The torquer was turned off after two 
periods (i.e., 4 pulses), and the state at that time was used <9 the initial 
state in order to compute the C-R bounds. One computer run was made for 
each (x, y and z) axis. The results revealed that the C-R bounds for modes 
with medium contribution to a particular axis gyro output (e.g. , mode 1 for 
x-axis in Table IV) are very large. For example, the C-R standard deviation 
at 130 sec. for using x-axis gyro data is 0.4 (p^ is O.Ol). Iherefore, 
only those modes which make strong contributions to a particular axis should 
be identified using gyro outputs along that axis. Figure (4) shows typical 
C-R standard deviations as they evolve with tine. Since the higher-frequen- 


1 

« 

i 
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cy node* decay much faster than the low-frequency nodes » the G-R standard 
deviations for the higher- frequency modes attain their steady»«state quickly* 
Frcn the data obtained, it can be concluded that the natural frequencies 
generally have the best identifiabililty (i.c., the smallest C-R standard 
deviations), followed by damping ratios, and then mode-slopes* Ihe 
excellent identif lability of the damping ratios is especially welcome, 
since it is usually the worst-known (apriori) of the parameters. Unlike the 
mode-slopes, the identifiability of the damping ratios (p.) improves as p.'s 

X 

decrease. 

The data generated is useful for comparing identifiability of different 
parameters rather than for obtaining the absolute identifiability. Trial 
application of algorithms such as the maximm likelihood method will reveal 
the actual accuracy with which the parsmeters can be identified. As men- 
tioned earlier, the G-R bounds for modes with weak or medium contributions 
to the sensor outputs are large. That is, those modes have poor identifi- 
ability with respect to sensor outputs along particular axes. Therefore, 
if one can compute the estimated accuracy of identification while perform- 
ing identification, one can conclude (assuming all modes are excited) that 
the mode-slopes for the poorly identifiable modes are small, and can be 
ignored (i.e., assumed to be zero) while designing the control systems. Ihe 
G-R bounds also indicate which sensor outputs shoul;ii be used to identify 
which modes. That is, the structure of the model used in identification can 
exclude modes which make weak or mediim contributions to the sensor outputs, 
and have poor identifiability. 

GONGLUDING REMARKS 

The identifiability of the modal parameters of large flexible 
structures was investigated. The expression for the optimal initial state 
which maximizes the trace of the Fisher information matrix (in the absence 
of persistent inputs) was obtained. Gramer-Rao (G-R) lower bounds were com- 
puted for the modal parameters of the 122 meter hoop/column antenna. It was 
found that modal frequencies generally have the lowest G-R bounds, followed 
by damping ratios and mode-slopes (or shapes) . 
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Lower solar panels 


Fi 


Figure 1. Hoop /Column antenna concept 


10 





OF S 

w POOR OUAU 


Evolution of C-R bounds (typical) 



C-R Std. daviation 
(<^2^) - 6 X lO"^ 

" 0.01 
(true value) 


0 50 60 70 8C 30 100 no 120 130 
TIME ( SEC) 

Examples of C-R bounds 


0.75 

.04 

.01 

0.65 

0.5E-2 

0.28 

D 

.006 

.01 

0.87 

0.37E-2 

0.57 

8.8 

.08 

.01 

1 . 1 

12 

0.24E-2 

15 




(cr% » C-R Standard deviation expressed as a percentage of the parameter) 


Figure 4. Cramer-Rao bounds 
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TABLE II. MODAL CONTRIBUTIONS TO OUTPUTS 
*: Strong /: Medium (blanlc) : Weak 





















